Nonlinear modes in finite-dimensional TT-symmetric systems 



(N 

o 

(N 
< 



PL, 

d 



> 

in 

m 

(N 
O 
(N 



D. A. Zczyulin and V. V. Konotop 

Centro de Fisica Teorica e Computacional and Departamento de Fisica, Faculdade de Ciencias, 
Universidade de Lisboa, Avenida Professor Gama Pinto 2, Lisboa 1649-003, Portugal 

By rearrangements of waveguide arrays with gain and losses one can simulate transformations 
among parity-time (PT-) symmetric systems not affecting their pure real linear spectra. Subject to 
such transformations, however, the nonlinear properties of the systems undergo significant changes. 
On an example of an array of four waveguides described by the discrete nonlinear Schrodingcr 
equation with dissipation and gain, we show that the equivalence of the underlying linear spectra 
implies similarity of neither structure nor stability of the nonlinear modes in the arrays. Even the 
existence of one-parametric families of nonlinear modes is not guaranteed by the VT symmetry 
of a newly obtained system. Neither the stability is directly related to the VT symmetry: stable 
nonlinear modes exist even when the spectrum of the linear array is not purely real. We use 
graph representation of PT-symmetric networks allowing for simple illustration of linearly equivalent 
networks and indicating on their possible experimental design. 

PACS numbers: 42.65.Wi, 63.20.Pw, 05.45.Yv, 64.60.aq 



The effect of dissipation or gain on dynamics of a phys- 
ical system is a fundamental issue either in classical or 
in quantum theories. Optics is one of the areas where 
the respective models appear naturally and are explored 
already for many years in contexts of different kinds of 
dissipative solitons One of a number of widely used, 
fundamental and simple models is an array of waveguides 
in the presence of gain and losses 0. This model is de- 
scribed by the discrete nonlinear Schrodingcr equation 
(DNLSE), which is fairly general. Its applications range 
from the so-called discrete optics Q to biophysics 4 1 and 
the meanfield theory of Bosc-Einstein condensate 5| (for 
a broad range of applications of DNLSE see also 6(). 

Recently, great interest in systems with dissipation 
and gain was triggered by the discovery of the so-called 
parity-time (PT) potentials, which in a definite range of 
parameters obey purely real spectrum Q . Numerous lin- 
ear physical systems for which VT symmetry is of great 
relevance have been proposed. Among them we men- 
tion non-Hermitian extension of quantum mechanics Q , 
electromagnetic wave propagation in a planar waveguide 
filled with active media and beam propagation in 
optical lattices [l(J. The phenomenon of VT symme- 
try breaking has been experimentally implemented in op- 
where the equations governing the system were 
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tics 

earlier known as describing a unidirectional coupler, i.e. 
as a particular form of the DNLSE . 

Nonlinear 'PT-symmetric problems were first posed in 
the context of the quantum field theory accounting for cu- 
bic interactions 12 [ and in guided wave theory |1 3| . Be- 
ing natural for optical applications, the nonlinear prob- 
lems received particular attention in the context of exis- 
tence of gap solitons 13 1 and defect modes Q in PT- 



symmetric lattices, as well as in context of the nonlinear 
PT-symmetric couplers in stationary 15, [l6| and soli- 
tonic regimes. More generally, the nonlinearity en- 
riches possible statements of the problem allowing for 



including the effects on nonlinear [18| . as well as both 
linear and nonlinear [l9j PT-symmetric potentials. 

It is known that by applying a similarity transfor- 
mation to a given linear PT-symmetric system, a new 
system with real spectrum can be constructed. Thus 
in [2(3] new potentials (not necessarily PT-symmetric) 
with real spectra were constructed using the Darboux 
transformation, while in plj pseudo-Hermitian opera- 
tors were introduced and unitary equivalence of PT- 
symmetric and Hcrmitian operators was established. It 
turns out, however, that possible mutual reductions of 
Hermitian, PT-symmetric, and pseudo-Hermitian linear 
operators leaving the spectrum pure real, may introduce 
dramatic changes in the properties of the respective non- 
linear systems. The analysis of such changes is the main 
goal of the present Letter. 

More specifically, we show that PT-symmetric systems 
obeying the same linear spectrum may either have one- 
parametric families of nonlinear modes or have not. If 
the families exist, stability of the modes is essentially dif- 
ferent for different systems, still having the same linear 
spectrum. Moreover, stable nonlinear modes may exist 
beyond the VT symmetry breaking. For a discrete sys- 
tem consisting of four waveguides we find that breaking 
of VT symmetry can occur in two different ways: the lin- 
ear spectrum acquires either two complex and two real 
eigenvalues, or all four eigenvalues become complex. Fi- 
nally, we represent each underlying linear system by a 
graph, allowing one to catalog different linearly equiva- 
lent PT-symmetric systems. 

We consider an array of N waveguides (sites) and de- 
note the field in the nth waveguide by q n (z), where z 
is the propagation distance. If each waveguide have dis- 
sipation or gain described by 7„, positive or negative, 
respectively, then the field propagation is governed by 
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the DNLSE 



N 



- K nm q m - \q n \ 2 q n - ijnQn, (1) 

m—l 



where q n = dq(z)/dz. Here we admit the existence of 
non-local coupling among the waveguides, which is de- 
scribed by the coefficients K nm = K mn = K* m which 
will be treated as entries of the real symmetric matrix 
K: K = K^, where is the Hcrmitian conjugate ma- 
trix. It is convenient to introduce diagonal matrices G = 
diag(7i,...,7iv) and F(q) = diagflgi | 2 , \q 2 \ 2 , ■ ■ ■ , \qn\ 2 ), 
which describe the dissipation and the nonlinear part of 
the system, respectively. Then the system ([T]) can be 
rewritten in the form 



iq 



-[H + F(q)]q, H = K + iG. 



(2) 



We search stationary nonlinear modes in the form q(z) 
w, where b is the propagation constant, and w 



w, W11C1C u 

(lUi, W2, ■ ■ • , wn) t solves the stationary DNLSE 
bw = [H + F(w)]w. 



(3) 



Requiring the spectrum of the linear problem 6w = Hw 
to be real, which is necessary for all linear modes to be 
propagating, we impose the constraint Xm=i 7n = 0. 

The matrix H is PT-symmctric if it commutes with a 
VT operator: [PT, H] = 0. Hereafter P is an orthogonal 
symmetric (and therefore Hcrmitian) matrix, and T is 
clement-wise complex conjugation: 7~q = q*. Using that 
= THT = PUP, we observe that the linear system 
iq = — Hq admits an integral of motion (see also (22[) 
Q = jt {Pq, q) , where the inner product is defined as 

Now we can specify the problem at hand: we consider 
the existence and stability of nonlinear modes of PT- 
symmetric lattices whose linear parts are related to each 
other by similarity transformations, all having the non- 
linearity of the on-site type. Such a statement is nat- 
ural for arrays of optical waveguides, since linear links 
among them can be arranged by assembling waveguides 
in different geometries, while the dissipation/gain and 
the nonlincarity are the characteristics of each partic- 
ular waveguide, which can be routinely controlled (see 
also Fig. [3]). One of our main findings is that linearly 
equivalent PT- symmetric lattices result in qualitatively 
different properties of their nonlinear extensions. 

VT -symmetric "quadrimer"- Since the gain and dis- 
sipation must compensate each other, the simplest mod- 
els allowing for nontrivial distribution of the dissipation 
have three or four waveguides. Below we concentrate on 
a quadrimer, respectively setting N = 4. We start by re- 
visiting recently considered in [16l | system with the next- 



neighbor interactions: K n 



1 . The corrcspond- 




FIG. 1: "Phase diagram" (PD) for the linear quadrimer Ho. 
The dark-grey diamond-shaped domain corresponds to un- 
broken VT symmetry; in the light-gray domains there are 
two real and two complex eigenvalues. In the white domains 
all eigenvalues are complex. In panels (a)-(e), correspond- 
ing to the points (a)-(e) in the panel (PD), we show families 
of nonlinear modes for: (a) 71,2 = 0.25; (b) 71 « —0.37, 
72 w 1.49; (c) 71 « 1.49, 72 » 0.36; (d) 71,2 = (e) 
71 = 0.1, 72 = 0.95. Stable (unstable) modes are shown by 
solid blue (dashed red) lines. 



with respect to Vo = [ ^ ^ ) (hereafter o\ 2 3 are the 

Pauli matrices and is the 2x2 zero matrix) . Depending 
on particular values of 71,2, three different situations are 
possible: (i) unbroken or exact VT symmetry, when all 
the eigenvalues b n , n = 1, 4, of H arc real; (ii) broken 
VT symmetry with two real and two complex conjugated 
eigenvalues (notice that this is possible only if 71 ^ 72); 
(iii) broken VT symmetry with all b n complex. Thus, 
the "phase space" (71,72) can be divided into three do- 
mains as it is shown in the phase diagram (PD) of Fig. [T] 
A feature of the phase diagram is the existence of the 
triple points Tj, j = 1, . . . , 4, where the three domains 
touch. The triple points correspond to values 71,2 for 
which b n — for each n = 1, ...,4. Depending on how 
7i ; 2 change in vicinity of Tj, either the "PT-symmetric 
phase or one of the VT symmetry broken phases arise. 

If Ho is exactly "PT-symmetric, then its linear eigen- 
states w are simultaneously the eigenstates of the corre- 
sponding VT operator, i.e. VqTw = w (up to irrelevant 
phase shift). It is natural to look for nonlinear modes 
that possess the same property: VqTw = w. Therefore 
we require w\ = w\ , W2 = , which reduces Eq. ([3]) to 



ing matrix, which we denote as Hq, is PT-symmetric 
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FIG. 2: Families of nonlinear modes of the Hermitian 
quadrimers, whose linear parts are described by Uh, chosen 
to have the same linear eigenvalues b n as the corresponding 
PT-symmetric quadrimers Ho illustrated in the panel (a), 
(b), and (c) of Fig. Q] 



We represent wi = Wie 1 * , where W\ and are 
real. Then Eq. (|4ap gives w 2 = W 2 e I$ , where 
W 2 — W\(b — W 2 — 171) is complex, and from 
Eq. gbj we obtain e" 2 ^ = /(Wi), where /(Wi) = 
(&W 2 - Wi - |W 2 | 2 W 2 - i72W 2 )/W 2 *. If a root of the 
equation |/(Wi)| 2 = 1 is found, then w\ and iu 2 can be 
readily obtained. Thus nonlinear modes of the quadrimcr 
correspond to the roots of a single equation |/(Wl)| 2 = 1 
with respect to one real unknown W\ . It is a purely tech- 
nical matter to reduce the latter equation to: Ps(W 2 ) = 
0, where -Ps(0 is an eighth-degree polynomial with real 
coefficients. Each positive root of Ps(£) corresponds to a 
nonlinear mode of the quadrimer. Since the roots depend 
continuously on b, the nonlinear modes constitute contin- 
uous families for fixed parameters of the system pij ] . As 
it is customary, such families can be represented on the 
plane (U, b) where U = \ Y^ti=i \ w n\ 2 is the total energy 
flow in the array. Panels (a)-(e) of Fig. Q] illustrate typi- 
cal examples of the families, as well as linear stability of 
the modes. When 7i j2 belong to the domain of unbroken 
VT symmetry [see Fig. Q](a)], one observes four families 
branching off from the linear limit, i.e. from the points 
b = b n , U = 0. In Fig. [T](a) there also exist families that 
can not be continued from the linear limit. In panels (b) 
and (c) we also address the points that belong to the do- 
main of unbroken VT symmetry but are situated closely 
to the triple points T lj2 . In these panels one observes that 
after the bifurcation from the linear limit, all four families 
rapidly lose stability and two of them cease to exist if U 
is sufficiently large. Comparing panel (a) with panels (b) 
and (c), we conclude that increase of 71,2, i.e. approach- 
ing the VT symmetry breaking boundary, is unfavorable 
for existence and stability of the modes. However, the 
most surprising fact, is that stable nonlinear modes can 
be found in the domains of broken VT symmetry. Both 
in panel (d), which addresses the case when the spectrum 
consists of two real and two complex eigenvalues, and in 
panel (e), i.e. when all the eigenvalues are complex, one 
can find stable modes. 

Hermitian quadrimer- If Ho is exactly VT- 
symmctric, then there exists a unitary matrix R, 
which transforms Hq to a Hermitian matrix Uh 12 il l : 



limit the modes in the array with gain and losses 
described by H have the same propagation constants 
as the modes in the array without gain and losses, 
which is described by Uh- Hence, for any 7 lj2 lying 
in the domain of unbroken PT-symmetry of Ho, one 
can introduce a new DNLSE iq = — [H# + F(q)]q 



[c.f. @]. Following 2l|, one can find H# explicitly 
and observe that all its nonzero elements are real and 
given by H ffj i 2 = H H ,2i = hi, H H ,u = Hhai = h 2 , 
Hij )2 3 = Hh.32 = with hj being dependent on 
71.2 ■ By construction, the matrix Hjj has the same 
eigenvalues as H for the given 7^2. 

Unlike in the PT-symmetric case, the modes of the 
nonlinear system with linear part described by H# can 
be searched as real-valued and either even or odd, i.e. 



solving the system bwi = /iiw 2 ± h 2 Wi + wf, bu>2 — 
hiWi ±h 3 w 2 +W2, where "+" ("— ") stays for even (odd) 
modes. This system is equivalent to a fourth-degree poly- 
nomial equation with respect to w\ . Families of even and 
odd nonlinear modes of the Hermitian quadrimer are il- 
lustrated in Fig. [2] Comparing Fig. Q] and Fig. [2J we 
observe that even if the matrices Ho and Uh have the 
same eigenvalues, the respective nonlinear systems show 
considerable differences in the properties of modes. The 
most visible differences are: (i) for the Hermitian system, 
the families bifurcating from the linear limit never close 
forming a saddle-node bifurcation [c.f. panels (b) and 
(c) in Fig. Q]and Fig. |5]; (ii) the leftmost family of the 
Hermitian system is always stable; (iii) in general, stable 
nonlinear modes of Eq. © with Ho and H# correspond 
to different values of the propagation constant b. 

"Generalized" quadrimer- Being of the dissipative 
nature, the considered above PT-symmetric quadrimer 
with linear part described by Ho possesses a property, 
usually typical for conservative systems - for the given 
parameters of the system (inter-site interactions K and 
dissipation 7i, 2 ) its nonlinear modes constitute contin- 
uous families rather than appear as isolated attractors. 
This peculiarity of nonlinear PT-symmetric systems was 
reported in several studies [l3, 2 18, 19| . Here we ar- 



gue that existence of the continuous families of nonlinear 
modes is not a typical property of PT-symmetric sys- 
tems. Specifically, the nonlinear PT-symmetric systems 
that admit the families of the modes appear as "isolated 
points" in a continuous set of generic PT-symmetric sys- 
tems. 

To this end we focus on the particular case 71 = 7 2 = 
7, i.e. G = diag(7, 7, — 7, — 7), and introduce an one- 
parametric family of matrices Hq(/3) = Kq(/3) + iG with 



K0G8) = 



1 

1 cos j3 

cos P sin 2/3 

\ -sin/3 cos2/3 



\ 

— sin /3 

cos 2/3 
-sin 2^ J 



H 



H 



H 



This means that in the linear and real parameter /3. One can ensure that Hq(/3) is 



4 



PT-symmetric with respect to Po(/3) = 



PCS) 
p(/3) 

where p(f3) = cos/3 ci + sin/3 0-3. For (3 = the ma- 
trix Ko(0) includes only the next- neighbor interactions, 
i.e. Ho(0) is merely the linear part of the PT-symmetric 
quadrimer studied above (with 71 = 72). Definition of 
Ho(/3) guarantees that its eigenvalues do not depend on 
(3. But the eigenvectors of Ho(/3) do depend on (3. 
H 



Next, using M = 



fi 



/here fj, 



o~z + ia 2 , one 



can generate a new matrix Hi(/3) = MHq(/3)M 1 , where 
HiG8) = KiG8) + iG, 



K1G8) = 



/ 1 -k- k + \ 

-1 fc_ -k+ 

—k- k- cos 2/3 sin 2/3 

\ fc+ -fc+ sin 2/3 -cos 2/3 / 



and fc± = ^sin(/3 ± f ). Then Hi(/3) is PT-symmetric 
with respect to Pi(/3) = MP (/3)M- 1 . Notice that the 
transformation M does not affect the dissipative compo- 
nent iG, which is the same both for Ho(/3) and Hi(/3). 
Obviously, the eigenvalues of the matrix Hi (/3) are the 
same as for Ho(/3) and also do not depend on (3. 

To give better physical insight into the systems 
Ho,i(/3), in Fig.[3]we introduce their weighted graph rep- 
resentation. The vertexes of graphs correspond to the 
sites q n while the edges (lines) represent inter-site cou- 
pling having weights equal to the values of the respective 
matrix elements: e.g. a line between the vertexes qi and 
qi corresponds to the elements K\ i = i^2,i of the matrix 
K. Each vertex is supplied by the sign "+" or "— " cor- 
responding to gain and dissipation. Wc notice that the 
loop edges, which describe the on-site interactions K„ tn , 
are not shown as being not relevant for the present con- 
sideration. The graph representation can be viewed also 
as indication on how one could place and connect the 
waveguides in an experiment in order to obtain the de- 
sirable PT-symmeric quadrimer. It is worth noting that, 
say the bottom graph in left column [i.e. Ho(7r/2)] can 
be reshaped into the line distribution of the waveguides 
similar to the graph Ho(0). 

Existence of nonlinear modes- Turning to nonlinear 
properties of the arrays, whose linear links are described 
by Ho.i(/3), let us suppose that the nth eigenstate of 
the underlying linear problem few = Ho.i(/3)w gives 
rise to a family of nonlinear modes. Then in the vicin- 
ity of the bifurcation point the nonlinear modes can be 
described using the expansion w = ew„ + o(e), and 
b = b n + e 2 b^ + o(e 2 ), where e is a small parameter, b n 

and w„ are the eigenvalue and the corresponding eigen- 

(2) 

vector of Ho(/3) [or Hi (/?)]. The coefficient fe„ can be 
readily found: b^ = (F(w„)w„, w* }/(w n , w*). This 
means that the bifurcation of nonlinear modes is possi- 
ble only if Imfe„ = for all n (we may conjecture that 
this condition is also sufficient for existence of the modes, 




P = tt/4 



/? = 7T/2 



FIG. 3: Graph representation of the systems Ho i i(/3) for 
generic and particular values of /3. In terms of the optical ap- 
plications, the circles with "+" [or "— "] represent waveguides 
with gain [or losses], while the lines indicate the directions, 
along which the coupling of the field has to be arranged. 



what was observed in all our numerical simulations). The 

(2) 

coefficient 6„ is easily computable. In Fig. 2] (left panel) 

Im&i 2) is plotted for H (/3). Only at /? = (3 k = nk/2 
(2) 

the coefficient b n becomes real for all n and the system 
Ho(/3) admits continuous families of nonlinear modes, 
while for all other (3 nonlinear modes bifurcating from 
the linear limit were not found. 

To understand peculiarity of the values /3k we notice 
that relation PTw„ = w„ ensures that the denominator 
in the formula for 6„ is real for any (3: (w„,w*) = 
(PTw„,Tw„) = (Tw„,PTw n ) = (w*,w„). Mean- 
while (F(w„)w n , w* } can have nonzero imaginary part. 

(2) 

Then the reality of the coefficient b,, is ensured by an 
additional constraint PT(F(w„)w„) = F(w„)w„, which 
is satisfied only for /3 = f3 k - 

For the system Hi (/3) the situation is similar - the fam- 
ilies of nonlinear modes exist only for /3 = irk/2 where 
Imbi? = for all n. In the right panel of Fig. [4] we 
show families of nonlinear modes of the array whose lin- 
ear part is described by Hi(0) with 7 = 0.25. Comparing 
the latter panel with Fig. [T] (a) (which also corresponds 
to 71,2 = 0.25), we again notice that, whereas the corre- 
sponding arrays have the same eigenvalues in the linear 
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FIG. 4: ImbZ> (n = 1, ...,4) vs 13 for H (/3) (left panel) and 
the families of modes of Hi(0) for 7 = 0.25 (right panel). 



limit, nonlinear modes of those arrays have essentially 
different properties. 

To conclude, we have considered nonlinear properties 
of different "PT-symmetric lattices (discrete nonlinear 
Schrodinger equations with gain and dissipation), whose 
linear parts are related by similarity transformations pre- 
serving the spectrum. Such systems describe, in par- 
ticular, arrays of optical waveguides with either gain or 
losses, which arc properly arranged in the space. Alterna- 
tively, a physical realization of the described phenomenon 
is possible in arrays of Bose-Einstein condensates loaded 
in multi-well potentials, provided the atoms are elimi- 
nated from given wells and are condensed in the other 
wells, simulating in this way losses and gain. 

On the case example of a PT-symmetric quadrimcr 
we have shown that the spectral equivalence of the un- 
derlying linear systems implies neither similarity of the 
nonlinear modes nor their stability properties. We have 
found that the existence of one-parametric families of 
nonlinear modes is not guaranteed by the VT symmetry, 
and appears as a peculiarity of a system rather than a 
general property. It was also found that the stability of 
nonlinear modes is not directly related to the VT sym- 
metry: stable nonlinear modes exist beyond the VT sym- 
metry breaking threshold. If the system includes two dif- 
ferent dissipativc coefficients, then the "phase diagram" 
of the PT-symmetric quadrimcr allows for existence of 
"triple" points, where three different phases meet. Fi- 
nally, we have shown that use of graph representation 
of "PT-symmetric networks gives straightforward indica- 
tion on their possible experimental design in optics, and 
provides graphical illustration of linearly equivalent net- 
works. 
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